A possible route to grow a (Mn:Si(i_ x )Ge :r )-based diluted magnetic semiconductor 
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We systematically investigate, using ab initio density-functional theory calculations, the properties 
of interstitial and substitutional Mn in both Si and Ge, as well as in the Sii-xGe^ alloy. We show 
that volume effects are not the main reason Mn prefers to be a subsitutional impurity in pure Ge, 
and chemical effects, therefore, play an important role. Using realistic models of Sii-zCe^, we show 
that for x > 0.16 substitutional Mn in Ge-rich neighborhoods become more stable than interstitial 
Mn, which may allow the growth of Si-based diluted magnetic semiconductors. 
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Diluted magnetic semiconductor (DMS) materials [l|, 
0, El- [J] have been know for a long time JaLhowever, it 
was the relatively recent growth of III-V [3, 111 Mn-based 
ferromagnetic DMS that have brought these materials 
to the spotlight, motivated by the possibility of fabricat- 
ing useful spintronic devices. Although these accomplish- 
ments are rather significant, if one considers the possible 
technological impact there is no doubt that the synthesis 
of a Si-based DMS material would be of great relevance. 
Ferromagnetism in the Mn^Gei-a; compound has been 
reported by more than one group 0, which makes one 
ask the obvious question: would it be possible to grow a 
similar Mn^Sii-a; ferromagnetic material? 

Considering this question in a recent paper ||, we 
have argued that: (i) as Si and Ge have similar elec- 
tronic structures, it is reasonable to expect that for large 
enough x, Mn^Sii-^ would become ferromagnetic below 
some critical temperature, in the same way Mn^Gei-^ 
does it. This hypothesis is supported by a recent cal- 
culation [9| where both MnuGe^j. and Mn^Sii-j; were 
studied. The question then is if the required values of x 
can be achieved; (ii) in Si, a Mn impurity favors an inter- 
stitial site, whereas in Ge it prefers a substitutional site 
(see below); (hi) as a consequence, a Mn substitutional 
impurity in Ge cannot diffuse as easily as an interstitial 
Mn in Si 0, fill fl2T ] , allowing the introduction of a large 
enough number of impurities without their diffusion and 
subsequent clustering. Therefore, it seems that a cru- 
cial issue to obtain a Mn^Sli-a material is to find a way 
to have the Mn occupying substitutional sites. We have 
previously found || that there are substitutional sites at 
the Si(100) surface with formation energies comparable 
to interstitial sites, opening up the possibility of placing 
Mns; through some non-equlibrium growth process, akin 
to the low-temperature molecular beam epitaxy proce- 
dures that are crucial to obtain III-V based DMSs 6]. 

In the present paper we address the following ques- 
tions: (i) Is it simply a volume effect that would make 
Mn prefer a tetrahedral interstitial site in Si whereas it 



prefers a substitutional site in Ge? (ii) Will a Ge rich 
environment in the SiGe alloy also favor a substitutional 
Mn site? Regarding question (i), we find that the cause 
Mn is substitutional in Ge is not simply a volume effect, 
and therefore finding ways to increase the Si lattice pa- 
rameter would not suffice to stabilize the substitutional 
Mn impurity. However, since specific chemical interac- 
tions are the possible explanation for the differences be- 
tween Mn in Si and Ge, this implies that the answer 
to question (ii) above is yes. We actually find that the 
substitutional Mn site shoud be the predominant one for 
Ge concentrations in the Sii-aGe^ alloy for x > 0.16. 
These findings have the important consequence that a 
Mn doped SiGe alloy may present similar ferromagnetic 
properties as the Mn^Gei-^ material, opening up the 
possibility of Si-based spintronics. 

All our results are based on ab initio calculations, 
based on spin-polarized density functional theory within 
the generalized gradient approximation (GGA) \V.\ . We 
have used ultrasoft pseudopotentials |bj] and a plane 
wave expansion up to 230 eV, as implemented in the 
VASP code 0. We have used a FCC-based supercell 
containing 128 sites, and the L-points for the Brillouin 
zone sampling. In all calculations the positions of all 
atoms in the supercell were relaxed until all the forces 
components were smaller than 0.02 eV/A. 

We initially considered pure Si and Ge crystals. We 
placed a Mn atom both at an interstitial as well as at 
a substitutional site. The formation energy of a neutral 
interstitial impurity [lfij . Ej, is calculated as 

Ef = Edcf — Eb u ik — fiMn (1) 

where Ed e f is the total energy of the supercell with the 
defect, Etuik is the total energy of the supercell without 
the defect (either pure Si or pure Ge), and /ijn„ is the 
Mn chemical potential For neutral substitutional 

impurities |l6j | the formation energy (Ef) is given by 

Ef = (Edef + Hx) — Ebulk — HMn (2) 
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TABLE I: Formation energies (in eV) for interstitial and sub- 
stitutional Mn in bulk Si and Ge. For both materials, the 
results are reported for the Si (ast) as well as the Ge (aGe) 
lattice constants. 

Lattice Constant MnJ 1 Mn s , Mn^ e Mn Ge 
aji 2J5 2^8 3l l7T~ 

a G e 2.5 2.7 2.1 1.5 



where [ix is the chemical potential of either Si or Ge [lj| • 
In Table [I] we present the results for Mn atoms in Si 
and Ge at their equilibrium lattice constants (we find 
a S i = 5.445 A for Si and a Ge = 5.750 A for Ge). As ex- 
pected 113 j f° r silicon the Mn atoms prefer to be at in- 
terstitial sites, with an energy difference Ai? s_/ between 
the substitutional and interstitial sites of Ai? s_/ = 0.3 
eV. For Mn in Ge, on the other hand, we obtain that the 
Mn substitutional impurity has a lower formation energy 
by AE S ~ ! — —0.6 eV. In order to investigate how much 
of this difference between Si and Ge is a result of their 
lattice parameter difference, i.e., how much is due to a 
volume effect, we repeated the same calculations for Si 
(Ge) in the lattice parameter of Ge (Si) [l7j . The re- 
sults are also presented in Table [I] In silicon, no sig- 
nificant changes were observed, with a small reduction 
in the substitutional Mn formation energy, which led to 
AE S - J = 0.2 eV instead of 0.3 eV. For Mn in Ge, the 
interstitial site becomes highly unfavorable, with the for- 
mation energy changing by more than 1 eV. The forma- 
tion energy for the substitutional site, on the other hand, 
changed very little. This leads to a large increase in the 
energy difference between the substitutional and inter- 
stitial sites: AE S ~ ! = —1.7 eV. Therefore, even though 
increasing the Si lattice parameter has a small tendency 
towards favoring a Mn atom in a substitutional site, the 
difference in lattice parameter between Si and Ge cannot 
account, by itself, for the distinct behavior of Mn in these 
materials. Thus, specific chemical interactions between 
the Mn atom and the host atoms are crucial in determin- 
ing the distinct behavior in Si and Ge. This suggests that 
in the Sii-^Ge^ alloy the Mn atoms may prefer to be at a 
substitutional site with a Ge rich environment, instead of 
at an interstitial site. In order to confirm this possibility, 
we performed calculations for both interstitial as well as 
substitutional Mn in Sii-^Ge^, for x=0.25, 0.5 and 0.75. 
Since the Sii-^Ge^ is a substitutionally disordered alloy, 
the vicinity of a Mn impurity is not uniquely determined, 
and for either the substitutional or the interstitial Mn, 
there are five different types of sites, if only the nearest- 
neighborhood is considered, i.e., a Mn surrounded by a 
configuration of Si and Ge atoms that can be labeled as 
Si l, Ge 4 ~ lJ , for v varying from to 4. 

To simulate the Sii-^Ge^ disordered alloys, we used 
supercells with 128 atoms where the atoms were dis- 
tributed as Special Quasirandom Structures (SQS) (lif . 
Details of the preparation of the supercells were given 



elsewhere [19(, and a similar procedure has already been 
used to study vacancies in this alloy [13, • It should 
be mentioned that due to the SQS approach, the disor- 
der of the alloy is taken into account in an explicit way. 
As we have shown before 0] , the alloy lattice parameter 
follows very closely the Vegard's law, and we therefore 
use asi^^c^ = (1 - x)a S i + xa Ge . 

To study the interstitial Mn defects, we randomly se- 
lected five sites with distinct first neighborhoods. As 
these sites will also have different second, third, etc., 
neighborhoods, an averaging procedure should be per- 
formed. From previous studies of vacancies in the alloy 
[20I I2H , we estimate an overall variation in the forma- 
tion energies of the order of ± 0.1 eV due the different 
vicinities, and hence this averaging will not alter our con- 
clusions in any significative way. The formation energy 
Ej(u, x) of a neutral interstitial [l6j Mn in the Sii-^Ge^ 
alloy, in a given neighborhood Si !y Ge 4 ~ l/ , labeled by v, is 
given by 

Ej(v,x) = Ej(v,x) - E bulk (x) - fi Mn (x). (3) 

Here Ej(v,x) is the total energy of the SQS structure 
with the Mn interstitital, Ef )U ii i .(x) is the total energy of 
the SQS alloy without any defects, and ^Mn(x) is the 
Mn chemical potential in the Sii-zCe^ alloy, which we 
determine [2^ as 

^Mn(x) = (1 - x)fJ,MnSi + X^MnGe ~ HSi (1 _ x) Ge x - (4) 

To study substitutional Mn defects a similar procedure 
as described above was employed. Since an atom from 
the original SQS structure must now be removed, there 
is the extra possibility of having the Mn replacing either 
a Si or a Ge atom. The formation energy Ej (X, v, x) of 
a neutral ^| Mn substituting an atom X (A^S^Ge), 
in the Sii-^Ge^ alloy, with a neighborhood Si^Ge 4- ", 
labeled by v, is given by 

Ef(X,i/,x)=Es(X,v,x)- E bu i k (x) + nx(x) - HMn(x). 

(5) 

Here Es(X, v, x) is the total energy of the SQS structure 
with the Mn in place of an atom X (X=Si,Ge), which has 
a chemical potential Hx(x) in the Sii_ I Ge :c alloy |23| . 

All our results are presented in TablelTTland Fig. ^ F° r 
the substitutional Mn calculations, we always removed a 
Si atom to place the Mn impurity, except for x = 0.5 
where a Ge atom was also removed. These latter results 
are shown in parethesis in Table [H] As can be seen, 
the formation energies are always very similar, indicat- 
ing that our conclusions do not depend on Mn replac- 
ing either a Si or a Ge atom. For interstitial Mn, the 
differences between the formation energies for the local 
configurations Si°Ge 4 and Si 4 Ge° are always 0.4 eV, for 
all x, with the lowest values occuring in Si-rich vicinities 
(Si 4 Ge°). For substitutional Mn this formation-energy 
spread is also independent of x, and has a value of 0.3 



TABLE II: Formation energies (in eV) for the interstitial and 
substitutional Mn in Si(i_ x )Gea;, for x = 0.25,0.5 and 0.75, 
and for the distinct neighborhoods Si"Ge 4_ ". For the substi- 
tutional cases, Mn was always replacing a Si atom, except for 
x — 0.5 where Mn replacing Ge atoms were also considered 



(results in parenthesis). 
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cV. Ge-rich neighborhoods, however, have lower forma- 
tion energies. An important result is the fact that the 
lowest formation energy for a substitutional Mn is al- 
ways smaller than the lowest formation energy for an 
interstitial Mn. The difference between these two low- 
est formation energies is 0.1 eV for x = 0.25, 0.4 eV for 
x = 0.5, and 0.5 eV for x = 0.75. This indicates that 
already for the Sio.75Geo.25 alloy, there are local config- 
urations that make substitutional Mn the lowest energy 
structure. Another significant result is the fact that for 
x > 0.5, all substitutional configurations have formation 
energies that are smaller than the smallest formation en- 
ergy for interstitial Mn. The overall picture of our results 
can be appreciated in Fig. ^ Curve fittings to both the 
average values (continuous curves) of the formation ener- 
gies as well as to the lowest values (long-dashed curves) of 
the formation energies are also presented in Fig. ^ The 
average values curves cross at x ~ 0.14, whereas the low- 
est values curves cross at a; ~ 0.18. These results indicate 
that for x > 0.16 ± 0.02 there should be Ge-rich neigh- 
borhoods in Si^^Ce^ where substitutional Mn atoms 
becomes more stable than interstitial Mn. 

At thermodynamic equilibrium, and assuming that the 
defects are independent of each other, the population of 
Mn interstitial defects in Sii-aGe^ is given by 

N 1 ^) = N s P v {x)exp{-G I } {v)/k B T), (6) 

where N$ is the total number of sites in the lattice, P v {x) 
is the probability |2J] for the v configuration to exist in 
the x-concentration alloy, G^(y) is the Gibbs free forma- 
tion energy of interstitial Mn at configuration is, and T 
is the temperature. Note that NsP v {x) is the effective 
number of interstitial sites of type v. Under the same as- 
sumptions, the population of Mn substitutional defects 
in Sii-zGe^ are given by 

N s (y) = N s P"(x)f(X, x)exp(-G s f (X, v)/k B T), (7) 

where f(X,x) = (1 - x) if AT=Si and f(X,x) = x if 
A"=Gc, and Gj(X, v) is the Gibbs free formation energy 
of substitutional Mn at configuration with Mn replac- 
ing an atom X (AT=Si,Ge). Assuming that, for a given 
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FIG. 1: Formation energies for interstitial and substitutional 
Mn in Sifi-^jGex, for x = 0.25, 0.5, and 0.75, and for the dis- 
tinct neighborhoods Si"Ge 4_ ". The results are slightly dis- 
placed in x (x — x ± 0.01) for clarity. Error bars (±0.1 eV) 
are estimates of the influence of vicinities beyond the nearest- 
neighborhood. Data for pure Si and Ge are also shown. 
The continuous lines are fits to average formations energies, 
whereas the long-dashed curves are fits to lowest formation 
energies. 

TABLE III: Relative populations for the interstitial and sub- 
stitutional Mn in Sifi-^Gez, for x = 0.25, 0.5 and 0.75, and 
for the distinct neig hborhoods Si"Ge 4_ ". The results were 
obtained according to Eq. |H1 using the data from Table [H] 
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temperature and composition, the vibrational entropies 
of all defects are similar and independent on their vicini- 
ties, we can calculate the relative population of the Mn 
defects as 

d / T x P v {x)exp[-Ef{u)/k B T\ 
11 Y, v , D P v {x)exp[-Ef{v)/k B TY W 

where D = I, S represents both the interstitial and sub- 
stitutional defects. For substitutional defects, since the 
formation energies for Mn replacing either a Ge or a Si 
atom are quite similar, and given that f(X,x) = 1, 
the above expression for n s (v,x,T) is obtained after a 
summation over X is performed. 

In Table II I II we present the relative populations for 
the interstitial and substitutional Mn impurities in the 
Si(i_ x )Gea; alloy, for x — 0.25, 0.5 and 0.75, for two tem- 
peratures, T=300 K and T=600 K. In all cases the overall 
population of substitutional Mn is larger than the inter- 
stitial one. At 300 K, even for x = 0.25 we already have 
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approximately 78 % of Mn at substitutional sites. As 
the temperature is increased to 600 K, this percentage 
decreases to 56 %. The significance of this result can be 
appreciated by noting that in bulk Si, at T=300 K, ba- 
sically 100 % of Mn are at interstitial sites, whereas at 
T=600 K there are approximately 99.7 % of interstitial 
Mn. This means that by alloying with Ge one can re- 
vert the overall population of defects from interstitial to 
substitutional Mn. 

In conclusion, we have shown that volume effects can- 
not be the main reason Mn prefers to be a subsitutional 
impurity in pure Ge. Chemical effects, therefore, must 
clearly play an important role. Through a systematic 
study of interstitial and substitutional Mn atoms in re- 
alistic models of the Si(i_ x )Ge x alloy, we have shown 
that for x > 0.16 the substitutional Mn in Ge-rich neigh- 
borhoods becomes more stable than the interstitial Mn. 
By playing with the temperature and the concentration 
x, and maybe also using non-equilibrium growth condi- 
tions H, it should be possible to obtain Si(i_ x )Gea; alloys 
with a large enough concentration of substitutional Mn 
atoms. Considering that in all calculations in the alloy 
we obtained a Mn local moment with the same value as 
in the pure crystals, and given that recent studies 
have shown that both Mn^Sii^a; and Mn^Gei-a, should 
have similar magnetic properties, all these facts indicate 
that Mn:Si( 1 _ :E )Ge 2 : is potentially a magnetic material 
like Mn^Gei-a, 0, opening up in this way the road to- 
wards Si-based DMS. 
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